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Abstract 

Based on a recent manifestly covariant time-ordered approach to the rela- 
tivistic many-body problem, the quark propagator is defined by a nonlinear 
Dyson-Schwinger-type integral equation, with a one-gluon loop. The resulting 
energy-dependent quark mass is such that the propagator is singularity-free 
for real energies, thus ensuring confinement. The self-energy integral con- 
verges without regularization, due to the chiral limit of the quark mass itself. 
Moreover, the integral determines the low-energy limit of the quark-gluon 

coupling constant, for which a value of <7 2 /4-7r = 4.712 is found. 
PACS number: 12.38.Aw 
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Typeset using REVTpjX 



In the absence of a direct derivation of the dynamical mechanism of quark confinement 
from the QCD Lagrangian, a number of different model approaches have been considered in 
recent years to achieve confinement (For a recent extensive review of related questions 

and further references, see ||.) Conceptually arguably one of the most elegant mechanisms to 
this end is based on the idea that the quark propagator should not have a pole, a mechanism 



which was proposed already twenty years ago [|K||lI]] anc ^ which has been implemented in 
Refs. [|T]-[|,[F]|| . In this letter, we want to present a model of quark confinement along similar 
lines, i.e., we will derive a quark propagator without any singularities for real energies. 

Our model is based on the cluster-dynamical approach to the relativistic many-body 
problem proposed in Refs. |1^,|T^]. This scattering formulation is a time-ordered one which 



is manifestly covariant under an off-shell modification Ai of the Lorentz transformation C 
The transformation (e', q') = Ai(e, q) of an arbitrary four- vector (e, q) describing a cluster 
with three-momentum q and off-shell energy e is defined by relating the three-momenta 
q' and q in two different frames by Lorentz transformations involving only the respective 
on-shell energies u' and u>, i.e., (u/, q') = £(cu,q), where cu = co(m,q) = (m 2 + q 2 ) 1 / 2 , etc., 
with m being the mass. The relation between off-shell energies e' and e is then defined by 

e — uj{m, q) = e' — u(m, q') , (1) 

similar to Galilei transformations in Euclidean space. The modified transformations Ai, 
therefore, leave the differences between on- and off- shell energies invariant, and reduce to £ 
when going on-shell. Performing then all internal integrations with the appropriate covariant 
integral measure [cf. discussion after Eq. (Q) below], the time- ordered cluster formalism of 



Refs. 1 12,13] is manifestly covariant under Ai and hence Lorentz-covariant for arbitrary on- 



shell matrix elements. We note that, while Ai appears as a nonlinear transformation of four- 



vectors (e, q), it can be can be understood [13] as a covariant four-dimensional projection of 



a linear five- dimensional extension of the Lorentz transformation C, i.e., Minkowski space 
is a hypersurface of an underlying five-dimensional manifold whose extra, fifth, dimension 
is interpreted as the off-shell energy variable with transformation properties defined in (|l]), 
while the usual Minkowski-space energy component is mapped onto the mass shell. 

The simplest possible way of defining a time-ordered quark propagator in the approach of 
Refs. [T2] , [r3"| ] is given by the one-gluon-loop expansion depicted in Fig. 1. The nonlinearity- 
structure of this integral equation is similar to a Dyson-Schwinger equation [p4|3, albeit in 
a time-ordered framework. In the center-of- momentum system (CMS), the resulting dressed 
(time-ordered) quark propagator is given by 

V(m,0) 

t(E -m + iO) - 



E — m + iO 

T(m, 0) 



(E - M + iO)Z -£(£) 



(2) 



where the first identity is a definition of the quark mass m; M is the bare mass, Z the 
renormalization constant and S the self-energy given by the gluon-loop bubble of Fig. 1. 
V{m, q) = {4 + m)/2m is the covariant on-shell spin-1/2 projector, i.e., g M = (uj, q). (We 
recall in this context that in the present formalism |12],[13| all spin degrees of freedom are 



described with on-shell energies, as a consequence of the particular off-shell continuation 
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described above.) Equation (0) defines the quark mass as a dynamical quantity. As an 
invariant, it must be an implicit function of itself, m = m(E — m), where E — m is an 
invariant according to ([I]), hence t(E — m + iO) is manifestly covariant. 

In the CMS frame, the explicit energy dependence |15[] of m(E — m) — > m(E) is readily 



found from (|2|). Requiring in particular that the mass m(E) be a positive semidefmite 
quantity for all energies and using the result (to be shown below) that £(-E) > for all 
energies, it follows that the renormalization constant Z must be unity, Z — 1, and hence 

m(E) - M = £(£) . (3) 

If we suppose for the moment that m(E) is indeed such that the quark propagator (Q) is 
free of singularities for real energies, i.e., that E — m(E) is always negative, then we may 
introduce a new function mo(E) by writing T,(E) = —itiq(E)/(E — m); we then find 




m(E) - M = o + \ ( — - — +m 2 {E) . (4) 



For large negative energies, m(E) approaches the bare mass M from above if we assume 
that i71q(E)/\E — M\ goes to zero, i.e., M provides a lower bound for the dynamical mass 
m(E). In the following, we will put the bare mass M = 0, in other words, we require that 
the dynamical quark mass approach its chiral limit for large negative energies, 




m(E) = j + ^{^) 2 + <(E) . (5) 

Our results will verify that m(E) has indeed the energy behavior given by this equation. 

To this end, we take the quark-gluon vertex to be undressed and given simply by gj tJ- X a /2, 
where g is the coupling constant and A a and 7 M are the usual color and Dirac matrices, 
respectively. From ([3]), with M = 0, it then follows, using the rules of Refs. [|T2| , |T3[] for 
evaluating the self-energy bubble E(i£) and keeping in mind the invariance property (0), 
that 

4 Q 2 f +co , f ,q m(E-e-u + m) N , s 

m(E)=i--^—- de d 3 q— -. r '- — , 6 

V ; 3(27r) 4 i-oo J H 2qu(m,q) (E - e - u + iO) (e - q + iO) V ; 



where the variables of Fig. 1 were used. Recall here that in the present approach [jT2|JT~3 



three-momenta are calculated with on-shell energy parameters in the Lorentz transforma- 
tions, and d 3 qm/2qu) is the corresponding covariant integral measure, with normalization 
factors 2|q| = 2q and u/m due to the gluon and the quark, respectively. Note that all 
masses here are expressed in the CMS notation as explicit functions of the correspond- 
ing off-shell energy; hence, in the integrand the quark mass and the associated on-shell 
17 1 energy are coupled functions of each other and implicit functions of themselves, i.e., 



m = m(E — e — uj + m) and uo = uj(m(E — e — uj + m), q). The remaining ratio arises from 
the energy convolution of the time-ordered quark and gluon propagators whose detailed 
variable dependence follows from the invariance requirement (|]). iV is the spinor matrix 
element, 

N = u(0,s)^^±^N^ lv u(0,s) , (7) 
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where = (uj, — q) is the on-shell four-momentum of the dressed quark in the loop of Fig. 
1 and u, u are Dirac spinors of the quark external to the self-energy loop with third spin 
component s. N^ u is given by the numerator of the time-ordered gluon propagator, 

G^(e-q + iO) = . (8) 

We treat the gluon here as an undressed particle similar to a bare photon. The question then 
arises how to incorporate gluon confinement. This is, of course, a very complicated and as 
yet unsolved problem (for recent reviews, see P,|16|, and references therein). We deal with it 
here by what we consider the simplest possible confinement mechanism known already from 
QED: discarding transverse contributions, we restrict the gluons to longitudinal and scalar 
(or, temporal) ones, which are unobservable the same way longitudinal and scalar photons 
are unobservable directly. We can offer no deeper justification of this procedure other than 
that it works for the present purpose. Hence, in Feynman gauge one has 

N" v = ^(3 V - rfT? , (9) 

where rf = (1,0) and /3 M = (0, n), with n being the unit vector defined by the gluon 
three-momentum q = q n. One finds then that simply N = — 1. 

In evaluating fl6|), we assume now that there are no singularities of the integrand in the 
lower half of the complex e plane other than the one at e = q — iO (cf. discussion below). 
The residue of the energy integration then yields 

. . q 2 [°° , m(E — q — uj + m) 1 . . 

m {E) = -f-J dqq ) — -L—- -— — — - ,(10) 

Stt z Jo uj (m(E — q — uj + m), q) E — q — uj (m(E — q — uj + m) , q) 

where we have also carried out the trivial angle integrations and dropped the iO in the 
denominator of the quark propagator in anticipation of the fact that the mass m(E) will 
indeed turn out to be confining. This equation is the main result of the present letter; it 
provides a highly nonlinear self-consistency equation for the dynamical quark mass m(E). 
Note that, in view of Eq. fl5|), for large momenta q the integrand behaves like m^{— 2q) / q 2 , in 
other words, the integral exists. [The numerical results discussed below verify that mjj(— 2g) 
goes to a constant for q — > 00.] The convergence of the integral thus is achieved by the chiral 
limit itself, without any cutoff whatsoever. 

Note also that Eq. fllOD does not fix any energy scale: using an arbitrary scale A, 
the corresponding solution m\(E) may be rewritten as m\(E) = h.fi\{E/K), where the 
functional form of the dimensionless fii is the same as m\. Furthermore, fii can be rescaled 
with an arbitrary (positive, dimensionless) constant A according to /ii(e) = A (JL\/\{e/ A). This 
is just the scaling behavior of the renormalization group. The energy scale must be fixed 
by other considerations, e.g., by the value of m(E) at E = 0, where m(0) = mo(0), which 
should be roughly given by the phenomenological constituent mass of the quarks, around 
300 MeV, which in turn is about equal to Aqcd- 

A trivial solution of Eq. (|10"D is the chiral limit itself, m(E) = 0. Since the right-hand 
side of Eq. (|T^) is a simple one-dimensional integral, other solutions can be obtained by 
an iterative numerical procedure. In view of the high nonlinearity, however, this is not as 
easy as it may look if one would like to find the complete set of solutions ||18||. At present, 
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therefore, we have made no attempt to do so but simply solved Eq. (|10D by straightforward 
iteration. As a starting point, we chose m(E) to be of the form (|5]) with m as a constant, 
m = 1, which fixes the energy scale. Both converged solution and starting function are 
given in Fig. 2. As can be seen, the solution always stays above the line given by the 
pole condition m = E. This solution thus is indeed confining, approaching the line m = E 
only asymptotically, as it must if quarks are to be free asymptotically. Note also that the 
starting function is an excellent approximation of m(E) for negative energies and around 
E = 0. Since Eq. ( |T0D is a nonlinear eigenvalue condition, the solution of Fig. 2 is found to 
correspond to a fixed coupling-constant eigenvalue of 

£ = 4.712 , (11) 

a value numerically accurate to better than 0.1%. With a straightforward iteration of Eq. 
(|10|), no other eigenvalue could be found. However, experience with iterative solutions of 



other nonlinear equations does not rule out the existence of other eigenvalues. 

As the low-energy limit of the QCD coupling constant, the value in (|TT| ) is a bit on 
the low side. We emphasize, however, that it was obtained without taking into account 
that the QCD coupling constant should be a function of the gluon momentum q in (|T0|). 
Solving Eq. ( |10|) with a decreasing running coupling constant a s (q 2 ) with a reasonable 
range of parameters, one finds, first, that a s (0) = g 2 /An may increase by up to about an 
order of magnitude — which does cover the range of expected low-energy values — and, 
second, that the corresponding solutions m(E) approach the asymptotic-freedom limit of 
m(E) = E faster than for the case shown in Fig. 1. In view of the relative simplicity of 
the present model, we do not want to dwell on any details of this investigation concerning 
a running coupling constant; nevertheless, we find it encouraging that the trends go in the 
right direction. 

Our result verifies that the dynamical quark mass can indeed be written in the form of 
Eq. (||). Since everything is completely known in the defining equation (|T0"D, in principle it 
should be possible to find the analytic properties of m(z) for an arbitrary complex energy 
z. In practice, this is complicated very much by the high degree of nonlinearity. Equation 
(|) suggests that m(z) has square-root cuts in the complex plane with branch points given 
by the implicit condition z = 2m(z). Various approximations and numerical studies seem 
to indicate that this condition cannot be met in the integrand of Eq. (|6|) for complex values 
of e in the lower half of the complex plane, thus justifying the derivation of Eq. fllPl). We 
mention, however, that we do not consider our findings to be entirely conclusive in this 
respect. 

In summary, we have presented here the simplest-possible one-gluon-loop expansion for 
the quark propagator in the covariant time-ordered formulation of the relativistic many- 
body problem of Refs. [|l2|,|ltj. We find a confining solution where the chiral limit of the 
dynamical quark mass m(E) itself provides the necessary cutoff to make the self-energy 
integral convergent. 

The author would like to thank Professors D. Schiitte and H. Petry for interesting dis- 
cussions. This work was supported in part by the U.S. Department of Energy, Grant No. 
DE-FG05-86-ER40270. 
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FIGURES 



FIG. 1. Nonlinear one-gluon-loop integral equation of the time-ordered dressed quark propaga- 
tor (solid line with dot) in terms of bare quark propagator (solid line) and gluon loop (cork-screw 
line); time proceeds from right to left. The off-shell energies and three-momenta are given in the 
center-of-momentum system of the dressed quark; e and q are the loop-integration variables of Eq. 

FIG. 2. Iterative solution of Eq. ([H]) (solid line). Both energy E — > e and mass m — > fi are 
mapped here by the same transformation x — > 2arctan(a;/s)/7r from (— oo,+oo) to the interval 
(— 1,+1); in the graph, the parameter s is arbitrarily chosen as s = 5. The energy scale is fixed 
by requiring m(0) = 1 [i.e., fi(0) = 0.1257 after transformation]. The long-dashed line is given by 
Eq. (j5|) with mo(E) fixed at mo(E) = 1. The short-dashed straight diagonal line corresponds to 
the pole condition m = E. 
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